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On 31 March 1983, the University of Delaware's Center for 
Remote Sensing Initiated a study to evaluate the spatial, radio- 
metric and spectral pe "mance of the Landsat Thematic Mapper for 
coastal and estuarine . ! es. The investigation was supported by 
Contract NAS-5-27580 from the NASA Goddard Space Flight Center. 

Cur research was divided into three major subprojects. A 
summary of the results for each subproject is given below. Details 
of each subproject are described in the three attached papers. 

1. A Comparison of Landsat TM to MSS Imagery for Detecting 

Submerged Aquatic Vegetation In Chesapeake Bay. 

Landsat Thematic Mapper (TM) and Mul tispectral Scanner (MSS) 
imagery generated simultaneously over Guinea Marsh, Virginia, were 
assessed in the ability to detect submerged aquatic, bottom-adhering 
plant canopies (SAV), , An unsupervised clustering, algorithm, was 
applied to both image types and the resulting classifications 
compared to SAV distributions derived from color aerial photography. 
Class confidence and accuracy were first computed for all water 
areas and then only shallow areas where water depth was less than 
six feet. In both the TM and MSS imagery, masking water areas 
deeper than six feet resulted in greater classification accuracy 
at confidence levels greater than 50t. Both systems performed 
poorly in detecting SAV with crown cover densities less than 70 %. 
On the basis of the spectral resolution, radiometric sensitivity, 
and location of visible bands, TM imagery did not offer a signi- 
ficant advantage over MSS data for detecting SAV in lower Chesapeake 
Bay. However, because the TM imagery represents a higher spatial 
resolution, smaller SAV canopies may be detected than is possible 
with MSS data. Detailed results of this investigation are described 
in the attached paper (Ackleson and Klemas 1985). 

2. Remote Sensing of Submerged Aquatic Vegetation: A 
Radiative Transfer Approach. 


Radiative transfer theory was used to model upwelling radiance 
for an orbiting sensor viewing an estuarine environment. The 
environment was composed of a clear maritime atmosphere, an opti- 
cally shallow estuary cf either clear or turbid water, and one of 
three possible bottom types: v-getation, sand, or mud. Upwelling 
radiance was calculated for each case in TM bands 1, 2 and 3 and 
MSS bands 4 and 5 using data available in the literature. A 
spectral quality index was usee to evaluate the relative effective- 
ness of TM and MSS bands in detecting submerged vegetation. 
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The effectiveness of an orbiting sensor in discriminating 
spectrally between submerged features is a function of the inherent 
contrast between the submerged features and how strongly the 
bottom signal is attenuated by the water column. In optically 
shallow water, the inherent contrast is the controlling factor. 
Thus, the optimum sensor band is that which correlates with the 
greatest inherent contrast between the submerged features. In the 
optically deeper water, the optimum sensor band is that for which 
the water column is most transparent. 

In the clear ocean water, the optimum band for detecting 
vegetation on a sand or mud bottom is shown to change with the 
optical depth of the water. In the turbid San Vicente water, the 
optimum band in optically shallow water remains the optimum band 
in optically deeper water. 

Under certain conditions, the spectral quality index will 

decrease' to zero for some intermediate depth and then increase, for 
yet deeper depths. This could have a significant impact upon how 
changes in water reflectance should be interpreted and can only be 
predicted with detailed knowledge of the spectral variation of 
optical parameters across the sensor band. 

Radiometric calibration is shown to be a dominant factor in 
determining which TM or MSS band will be optimum for detecting a 
submerged feature such as SAV . Sensors with higher gain settings, 
in general, produce data with more information. The gain built 
into the TM sensors are significantly higher than those of the 
MSS, giving TM data a clear advantage for detecting submerged 
features. In several cases, a band is found to be optimum for 
detecting SAV in a radiometric sense, anu yet sub-optimal in a 
spectral context. 

The Thematic Mapper appears to be useful for detecting 
narrow SAV beds In the case of Broad Creek, Maryland. A nonpara- 
metric classification of a 2 November 1982 image clearly discrimi- 
nated between SAV and the surrounding unvegetated bottom. 

3. Remote Sensing of Coastal Wetland Biomass Using Thematic 

Mapper Wavebands. 

Spectral data, simulating Thematic Mapper bands 3, 4 and 5 
were gathered in salt and brackish marshes using a hand-held 
radiometer. Simple regression models were developed equating 
spectral radiance indices with total live biomass for al terni- 
flora in a salt marsh and for a variety of plant species in a 
brackish marsh. Models were then tested using an independent set 
of data and compared to harvest estimates of biomass. In the salt 
marsh, biomass estimates from spectral data were similar to har- 
vest biomass estimates during mo3t of the growing season. Estimates 
of annual n^t aerial primary productivity calculated from spectral 


data were within 2M of production estimated from harvest data. 
During August, biomass estimates from spectral data in the brackish 
marsh were s'milar to t ' < omass estimated by harvesting techniques. 
At other times during the growing season, spectral data estimates 
of biomass were not always comparable to harvest biomass estimates. 
Reasonable estimates of wetlands biomass are possible during the 
peak of the growing season (August) using spectral data similar to 
Thematic Mapper bands 3, and 5 gathered with hand-held radio- 
meters . 

fr. General Results 

Eleven journal articles and proceedings papers were produced 
during this research. They are listed on the following page. A 
detailed paper on each of the three subprojects is included in 
this report. During the course of this investigation, additional 
comparisons of TM and MSS were made for coastal applications. 
Results of a comparison of Landsat MSS, TM and simulated SPOT data 
for coastal application were presented at symposia and are described 
in several articles (Aokleson et al. J9B5). , Some of t ,he re.sults 
of this comparison are summarized in Table 1. Basically, the 
better spatial resolution of TM and SPOT offer major improvement 
for detecting many features in the coastal zone which tend to be 
narrow and long. This includes beds of submerged aquatic vegeta- 
tion, submerged sand bars, pollution plumes, plots of marsh vege- 
tation, etc. The Thematic Mapper offers additional improvements 
due to its superior spectral bands. TM band 1 is particularly 
Important to studies of water properties and submerged features, 
while TM band 7 is very sensitive to moisture content of vegetation, 
which can be used as an indicator of plant vigor or stress. We 
have also found that all major wetland vegetation species can be 
clearly discerned in TM imagery. The spatial resolution of TM 
data appears to be better than 30 meters, i.e. it seems to be 
closer to 2b meters than 30 meters. 


TABLE 1. 

COMPARISON OF LANDSAT MSS 
AND TM FOR COASTAL STUDIES 
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APPLICATION 

IMPROVEMENT 

Vegetation and Land Use 
Mapping 

Medium 

Biomass Measurement 

Major 

Submerged Aquatic Vegetation 

Major 

Eutrophication 
Phytoplankton Blooms 

Medium 

Suspended Sediment Currents, 
Turbidity Fronts 

Medium 

Pollution Plumes 

Ocean Dumping 

Medium 

Bathymetry 
Erosion Control 

Major 

Ship Traffic 
Harbor Planning 

Major 

Gross Coastal Geomorphology 

Minor 

Sargassum 

Open Ocean Fronts 

Medium 

Internal Waves 
Sea State 

Minor 
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Summary 

On 31 March 1983, the University of Delaware’s Center for 
Remote Sensing initiate* a study to evaluate the spatial, radio- 
metric and spectral per nance of the Landsat Thematic Mapper for 

coastal and estuarine s es. The investigation was supported by 
Contract NAS-5- 27580 from the NASA Goddard Space Flight Center. 

Our research was divided into three major subprojects. A 
summary of the results for each subproject is given below. Details 
of each subproject are described in the three attached papers. 

1. A Comparison of Landsat TM to MSS Imagery for Detecting 

Submerged Aquatic Vegetation In Chesapeake Bay. 

Landsat Thematic Mapper (TM) and Multispectral Scanner (MSS) 
imagery generated simultaneously over Guinea Marsh, Virginia, were 
assessed in the ability to detect submerged aquatic, bottom-adhering 
plant canopies (SAV), , An supervised clustering . algorithm, was 
applied to both image types and the resulting classifications 
compared to SAV distributions derived from color aerial photography. 
Class confidence and accuracy were first computed for all water 
areas and then only shallow areas where water depth was less than 
six feet. In both the TM and MSS imagery, masking water areas 
deeper' than six feet resulted in greater classification accuracy 
at confidence levels greater than 50$. Both systems performed 
poorly in detecting SAV with crown cover densities less than 70$. 
On the basis of the spectral resolution, radiometric sensitivity, 
and location of visible bands, TM imagery did not offer a signi- 
ficant advantage over MSS data for detecting SAV in lower Chesapeake 
Bay. However, because the TM imagery represents a higher spatial 
resolution, smaller SAV canopies may be detected than is possible 
with MSS data. Detailed results of this investigation are described 
in the attached paper (Ackleson and Klemas 1985)* 

2. Remote Sensing of Submerged Aquatic Vegetation: A 

Radiative Transfer Approach. 

Radiative transfer theory was used to model upwelllng radiance 
for an orbiting sensor viewing an estuarine environment. The 
environment was composed of a clear maritime atmosphere, an opti- 
cally shallow estuary of either clear or turbid water, and one of 
three possible bottom types: vegetation, sand, or mud. Upwelling 
radiance was calculated for each case in TM bands 1, 2 and 3 and 
MSS bands 4 and 5 using data available in the literature. A 
spectral quality index was used to evaluate the relative effective- 
ness of TM and MSS bands in detecting submerged vegetation. 
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The effectiveness of an orbiting sensor In discriminating 
spectrally between submerged features is a function of the inherent 
contrast between the submerged features and how strongly the 
bottom signal is attenuated by the water column. In optically 
shallow water, the inherent contrast is the controlling factor. 
Thus, the optimum sensor band is that which correlates with the 
greatest inherent contrast between the submerged features. In the 
optically deeper water, the optimum sensor band is that for which 
the water column is most transparent. 

In the clear ocean water, the optimum band for detecting 
vegetation on a sand or mud bottom is shown to change with the 
optical depth of the water. In the turbid San Vicente water, the 
optimum band in optically shallow water remains the optimum band 
in optically deeper water. 

Under certain conditions, the spectral quality index will 
decrease to zero for some intermediate .depth and then increase, for 
yet deeper depths. This could have a significant impact upon how' 
changes in water reflectance should be interpreted and can only be 
predicted with detailed knowledge of the spectral variation of 
optical parameters across the sensor band. 

Radiometric calibration is shown to be~a dominant factor in 
determining which TM or MSS band will be optimum for detecting a 
submerged feature such as SAV. Sensors with higher gain settings, 
in general, produce data with more information. The gain built 
into the TM sensors are significantly higher than those of the 
MSS, giving TM data a clear advantage for detecting submerged 
features. In several cases, a band is found to be optimum for 
detecting SAV in a radiometric sense, and yet sub-optimal in a 
spectral context. 

The Thematic Mapper appears to be useful for detecting 
narrow SAV beds in the case of Broad Creek, Maryland. A nonpara- 
metric classification of a 2 November 1982 image clearly discrimi- 
nated between SAV and the surrounding un vegetated bottom. 

3. Remote Sensing of Coastal Wetland Biomass Using Thematic 

Mapper Wavebands. 

Spectral data, simulating Thematic Mapper ban's 3, it and 5 
were gathered in salt and brackish marshes ut. ng a hand-held 
radiometer. Simple regression models were developed equating 
spectral radiance Indices with total live biomass for S^ alternl- 
flora in a salt marsh and for a variety of plant species in a 
brackish marsh. Models were then tested using an independent set 
of data end compared to harvest estimates of biomass. In the salt 
marsh, biomass estimates from spectral data were similar to har- 
vest biomass estimates during most of the growing season. Estimates 
of annual net aerial primary productivity calculated from spectral 


data were within elf of production estimated from harvest data. 
During August, bioeass estimates from spectral data in the brackish 
marsh were similar Lr uiuu.ass estimated by harvesting techniques. 
At other times during the growing season, spectral data estimates 
of biomass were not always comparable to harvest biomass estimates. 
Reasonable estimates of wetlands biomass are possible during the 
peak of the growing season (August) using spectral data similar to 
Thematic Mapper bands 3* ** and 5 gathered with hand-held radio- 
meters. 

A. General Results 

Eleven Journal articles and proceedings papers were produced 
during this research. They are listed on the following page. A 
detailed paper on each of the three subprojects is included in 
this report. During the course of this investigation, additional 
comparisons of TM and MSS were made for coastal applications. 
Results of a comparison of Landsat MSS, TM and simulated SPOT data 
for coastal application were presented at symposia and are described 
In several articles (Aokleson et al. . 1 985). Some of t,he results 
of this comparison are summarized In Table 1. Basically, the 

better spatial resolution of TM and SPOT offer major improvement 
for detecting many features in the coastal zone which tend to be 
narrow and long. This includes beds of submerged aquatic vegeta- 
tion, .submerged sand bars, pollution plumes, plots of marsh yege- 
tation, etc. The Thematic Mapper offers additional improvements 
due to its superior spectral bands. TM band 1 is particularly 
important to studies of water properties and submerged features, 
while TM band 7 is very sensitive to moisture content of vegetation, 
which can be used as an indicator of plant vigor or stress. We 
have also found that all major wetland vegetation species can be 
clearly discerned in TM imagery. The spatial resolution of TM 
data appears to be better than 30 meters, i.e. it seems to be 
closer to 25 meters than 30 meters. 


TABLE 1. 


COMPARISON OF LANDSAT MSS 
AND TM FOR COASTAL STUDIES 


APPLICATION 


IMPROVEMENT 


Vegetation and Land Use 
Mapping 


Medium 

Biomass Measurement 


Mayor 

Submerged Aquatic Vegetation 


Major 

Eutrophication 
Phytoplankton Blooms 

- 

Medium 

Suspended Sediment Currents, 
Turbidity Fronts 


Medi urn 

Pollution flumes 
Ocean Dumping 


Medium 

Bathymetry 
Erosion Control 


Major 

Ship Traffic 
Harbor Planning 


Major 

Gross Coastal Geomorphology 


Minor 

Sargassum 
Open Ocean Fronts 


Medium 

Internal Waves 
Sea State 


Minor 
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ABSTRACT 


Landsat Thematic Mapper (TM) and Multi spectral Scanner 
(MSS) imagery generated simultaneously over Guinea Marsh, Virginia 
are assessed in the ability to detect submerged aquatic, bottom- 
adhering plant canopies (SAV). An unsupervised clustering algorithm 
is applied to both image types and the resulting classifications 
compared co SAV distributions derived from color aerial photography. 
Class confidence and accuracy are first computed for all water 
areas and then only shallow areas where water depth is less than 6 
feet. In both the TM and MSS Imagery, masking water areas deeper 
than 6 Ft. resulted in greater classification accuracy at confidence 
levels greater than 50Y. Both systems perform poorly in detecting 
SAV with crown cover densities less than 70X. On the basis of the 
spectral resolution, radiometric sensitivity, and location of visible 
bands, TM imagery does not offer a significant advantage over MSS 
data for detecting SAV in Lower Chesapeake Bay. However, because 
the TM imagery represents a higher spatial resolution, smaller SAV 
canopies may be detected than is possible with MSS data. 
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INTRODUCTION 


Submerged aquatic vegetation (SAV) is believed to play 
a major role in the ecosystem of coastal, estuarine, and inland 
waters. In Chesapeake Bay, species such as Zostera marina (eel 
grass) provide food, shelter, and breeding areas for waterfowl, 
fish, shellfish, and many other forms of aquatic life. Because of 
the enormous commercial value of these areas, there exists a need 
to periodically assess the distribution and abundance of submerged 
aquatic plant communities. 

' - 1 * ‘ • * *\H ‘ , V : 
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Historically, data concerning SAV distribution and abundance 
has been acquired through exhaustive field sampling programs. Today, 
with the exception of very small scale studies, such surveys have 
become prohibitively expensive. 

More recently, color aerial photography has been shown to 

provided much useful information for mapping SAV (Benton and Newnam, 
1976; Austin and Adams, 1978; Macomber, 1981; Orth et al., 1979; 
Orth and Moore, 1981 and 1983). Orth et al. (1979) and Orth and 
Moore (1981) were able to document changes in the distribution and 
abundance of SAV at six sites within Chesapeake Bay from aerial 
photography dating back to 1937. Orth ( 1 98^) was able to map SAV 
density within Lower Chesapeake Bay from color aerial photography 
collected in 1979, 1980, 1982, and 198*1. For each of these years, 
SAV distribution has been delineated oh 1:2*1000 USGS topographic 
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quadrangles and classified according to apparent crown cover density; 
<10%, 10-40%, 40-70% and 70-100*. 

Most recently, remote sensing data gathered from orbiting 
platforms such as the LANDSAT series satellites has shown promise 

in delineating SAV. The degree to which these data may be used 
depends upon the configuration of the particular sensor; i.e. spatial 
and spectral resolution, band location, and radiometric sensitivity. 
Jensen et al. (1980) successfully used Landsat Multispectral Scanner 
(MSS) imagery to map kelp ( Macrocystis pyrifera ) along the California 
coast. The spatial resolution of MSS data is approximately 60 
meters. Although the Landsat MSS collects data in four bands, two 
visible bands and two near infrared bands, only band 1 (500-600 
nm) and band 2 (600-700 nm) are useful for detecting submerged 
features because water is an efficient absorber of IR radiation. 

In 1983. NASA launched Landsat A equipped not only with a 

MSS hut the first of a new series of scanners, the Thematic Mapper 
(TM). The spatial resolution of TM imagery is approximately 28 
meters and 3 of the 7 available bands are located in the optical 
portion of the spectrum; band 1 (^50-520 nm), band 2 (530-610 nm) 
and band 3 (620-690 nm). In addition, the radiometric sensitivity 
of the TM is as much as 3 times that of the MSS (Table 1). All of 
these attributes suggest that TM imagery should contain more sub- 
merged features information than MSS data. The purpose of this 
work is to compare the relative capability of TM imagery with MSS 
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data for detecting SAV within Lower Chesapeake Bay. The compari- 
son Is based '«nnr SAV distribution derived by Orth from color aerial 
photography collected 28 days prior to the LANDSAT overpass. 


DESCRIPTION OF STUDY AREA 

The area of interest is located in Lower Chesapeake Bay, 
at the mouth of the York River, along the southern shore of Mobjaek 
Bay (Figure 1). The area encompasses 41.1 km*, approximately 6.5 
km* of which is inhabited by SAV. These bottom-adhering canopies 
are composed mostly of Zostera marina and Ruppla marltlma and are 
typically situated in shallow water adjacent to the shoreline. 
During the summer, both of these species are abundant with Ruppla 
occupying the most shallow areas (<.5m) and Zostera dominating 
deeper areas to a depth of about 1.5m. During the winter, the SAV 
distribution and abundance decreases as the Ruppi a practically disap- 
pears and the Zostera dies back significantly. 

The bottom substrate throughout the study area is composed 
of bright sand. From the air, the sav canopies appear as dark 

patches upon an otherwise bright bottom. 

The water clarity varies seasonally with most clear condi- 
tions occurring in the winter months. During the summer, higher 
nutrient concentrations, warmer water, and increased surface light 
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Intensities, often Induce plankton blooms. The resulting turbidity 
may greatly reduce the ability of an observer or remote sensing 
device to detect SAV. Van Tine and Wetzel (1983) reported seasonal 
averages of lrradlance attenuation measured at several sites In 
Lower Chesapeake Bay containing SAV. Their data indicate that 
during the winter, the average depth at which the downwelling light 
Intensity has decreased to 1J of the surface intensity for blue 
(1110 nm), green (5^0 nm) and red (671 nm) light Is, respectively, 
2.0 m, 5.7 m, and 3.2 m. During the summer, the mean attenuation 
depth decreases to about 0.8 m, 1.9 m„ and 1.6 m for blue, green 
and red light, respectively. 

The tidal range throughout the Guinea Marsh area Is approx-, 
imately 1 meter. Given the relatively turbid water conditions, 
especially during the simmer, remote sensing of SAV should be con- 
ducted as close to low tide as possible. 


METHODOLOGY 


Landsat TM and MSS imagery generated simultaneously over 
the study area was classified using an unsupervised clustering algor- 
ithm. All data analysis was performed using an ERDAS H00 computer 
system. The input parameters governing the operation of the clus- 
tering program are defined as functions of the TM and MSS system 
noise. The classified TM and MSS imagery are then compared, pixel- 


by- pixel, to digitized SAV distributions delineated from color 
aerial photography. Class confidence and accuracy are calculated 
for eacn scene and used as a basis for comparing the two image types. 
LANDSAT Imagery 


LANDSAT TM and MSS Imagery of the Guinea Marsh area was 
generated simultaneously on July 19, 198 1 *. At the time of the 
overpass, 1009 local time, the sky was cloud-free and reasonably 
clear of haze. The solar zenith angle was 32 s . The tide at Guinea 
Marsh was flooding and low tide occurred at 0706 local time. 

Prior to classification, each image was geometrically cor- 
rected and registered to a USGS topographic quadrangle of the Guinea 
Marsh area. The TM image was resampled to simulate 28 meter reso- 
lution and the MSS data was resampled to simulate 60 meter resolu- 
tion. 


Reference SAV Maps 


Since 1978, the Virginia Institute of Marine Science has 
conducted annual and biannual surveys of SAV distribution within 

Lower Chesapeake Bay (Orth and Moore, 1983). In each sur vey, 1:2*1000 
color aerial photography was collected during the summer months. 
The photography was then classified into several catagorles of SAV 
using visual interpretation techniques. SAV catagorles are based 
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upon spatially averaged crown densities and expressed as a percent 
coverage or crown cover; 0-10%, 10-40%, 40-70%, and 70-100%. The 
SAV delineations were then transferred to 1:24000 USGS topographic 
quadrangles. The 1984 survey, conducted on June 21, 1984, 28 

days prior to the Landsat Overpass, is used here as a reference 
with which to compare the classified TM and MSS imagery. A TM 
reference image is created by digitizing the SAV map, using a CALCOMP 
9100 digitizer, to simulate 28 meter resolution. To create the MSS 
reference image, the SAV map is digitized to simulate 60 meter 
resolution. In both cases, the lmage/SAV map registration accuracy 
is it pixel. . 

Landsat TM And MSS Image Clustering 

The image clustering algorithm, CLUSTER, was developed by 
ERDAS Inc. and is based upon a program formulated by NASA/JSC as 
part of the ASTEP software package ( ERDAS 400 Image Processing and 
Geographic Information System User's Guide , 19Cj1» CLUSTER is a 
two-pass unsupervised classifier designed to group pixels of a 
multi-band image into N number of spectral classes. In the first 
pass, the operator defines the maximum cluster radius (R), the 
minimum Euclidean distance in P-space between cluster centroids 
(D), where P is the number of spectral bands in the image, and the 
maximum number of classes to be produced (N). The program considers 
each pixel in a sequential manner, starting at the upper left corner 
of the scene and progressing line-by-line to the lower right corner. 
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Class mean vectors, generated in the first pass, are used In the 
second pass, wherein a minimum distance classifier is applied to 
the entire image. 

One disadvantage to using an unsupervised classifier such 
as CLUSTER is the relatively small amount of control that the oper- 
ator has over the program. Because image classification is a f'inc- 
tion of the input variables, the operator must have some prior 
knowledge of the scene variance in order to select reasonable values 
of R, D and N. In this work, it is also necessary to standardize 
the selection of input values in order to provide a basis for com- 
parison between the classified TM and MSS imagery. 

As an image is generated by either the TM or the MSS, system 
noise is introduced, primarily as a result of differences in detector 
calibration, voltage fluctuations within the satellite, and round 
off error in the digitization process. The effect of this noise 
is to limit the ability of the sensor to detect more subtle varia- 
tions in scene radiance than would normally be detected if the image 
were noise-free. Within any classification, no matter how sophis- 
ticated the classifier, classes will be meaningful only if the 
spectral separation from other classes is greater than the system 
noise. 


The maximum cluster radius and minimum distance between 
dusters are defined as functions of the system noise as 




-rr 


MS A.-*,.. 3. 
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R ’ ci s ( 1 ) 

and 

D - 2R , (2) 

where o 3 a is the variance in count values as a result of system 
noise. The system noise variance may be estimated for each image 
band by calculating o* within a homogeneous region of the image. 
The more homogeneous the region and the larger the number of pixels 
within the region, the better the estimate of o s 2 . Within both 
the TM and MSS imagery, optically deep water appeared to be most 
homogeneous and geographically similar areas were selected within 
each image for the estimate of o 3 l . Within the TM image the area 
comprised 90 pixels while, in the MSS image, 20 pixels were repre- 
sented. Values of R and D used to classify the TM and MSS imagery 
are shown in Table 2. 

For each classification, N-50, the maximum number of classes 
that CLUSTER can create. However, classification error can occur 
when actual number of classes created is N. This will most likely 
be the case when the selected values of R and D are small relative 
to the variance in count values throughout the scene. When the 
maximum number of clusters is reached, new clusters may not be 
created until two or more of the existing classes are aggregated. 
Within the TM and MSS scenes, the majority of scene variance occurs 
within land areas. Therefore, prior to running the classifier, 
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all land areas within both images were masked using one of the IR 

bands - TM band 9 '770-900 nm) and MSS band 9 (800-1100 nm). This 
caused the classifier to assign just 1 class to land and allowed 
as many as 99 water classes to be created. 

Calculation Of Confidence And Accuracy 


The classified TM ant MSS images are compared pixel-by- 
pixel with the digitized, 1989 SAV distribution map of Guinea Marsh. 
For each class, confidence and accuracy are calculated, where 


Confidence 


# Correctly Classified Pixels 
# Pixels Within Class 


(3) 


and 


# Correctly Classified Pixels 

Accuracy - # SAV Map Pixels 


(9) 


Within each class, a pixel is correctly classified if it correlates 

with a SAV pixel within the reference image. 


TM and MSS classification confidence and accuracy was first 
computed for all water areas. Next, all water areas deeper than 6 
feet were masked within the TM and MSS classifications using N0AA 
chart 12238, published in 1979, digitized and registered to each 
image. Class confidence and accuracy was then recalculated and 
compared with the values representing all water area3. 
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RESULTS AND DISCUSSION 

Classification accuracies computed for all water areas 
indicate no significant difference between TM and MSS imagery for 
detecting SAV in the Guinea Marsh area (Table 3)« Both data types 
performed poorly in discriminating SAV canopies having crown cover 
densities of 0-1 0$, **0-70%, and 70-1 00$. In the latter case, clas- 
sification accuracies appeared to increase within the 50 - 75 % and 
75-100$ confidence intervals, but the increases were too small to 
be considered significant. 

There appeared to be moderate success in detecting SAV crown 
cover densities of lO-^Ot. This was a somewhat surprising result 
as classification accuracies were expected to increase with SAV den- 
sity. To further investigate this result, all pixels having a 
classification confidence 250$ were highlighted within the unclas- 
sified TM and MSS imagery. A high correlation was found between 
these pixels and bright areas within the imagery. 

For any plant canopy, the volume reflectance may be divided 
into 2 components; that light which is reflected from the plant 
material and that light which is reflected from the underlying 
substrate. When the canopy is dense (i.e. crown cover densities 
a pp: oaching 100$), the first component dominates the volume reflec- 
tance of the canopy. As the plant density decreases, the second 
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component Increases. Within the Guinea Marsh area, the substrate 
is composed of bright sand and dense SAV canopies appear from the 
air as dark areas upon an otherwise bright bottom. The fact that 
pixels classified with high confidence for the 10-^0% SAV crown 
cover correlate with bright areas within the unclassified imagery 
indicates that the classifier is not discriminating SAV, but rather 
shallow sandy areas. Although aerial photography indicated a 1CM0J 
SAV crown cover, the TM and MSS scene radiance is dominated by 
reflectance from the bright, sandy substrate. 

When all water areas of depth *6 feet are eliminated from 

the classification analysis, accuracies for the 50-75% and 75-100% 
confidence intervals were found to increase only for SAV of crown 
cover density 70-100 % (Table 4). For the 75-100% confidence inter- 
val, TM classification accuracy increased from 0.2% to 19.7%. 
MSS classification accuracies increased similarly; 0.5% to 13-3%. 
All other accuracies remained unchanged for confidence levels greater 
than 50%. This indicates that within both TM and MSS imagery, it 
is difficult to distinguish between Lower Chesapeake Bay SAV and 
optically deep water. More Importantly, by masking optically deep 
water, classification accuracies associated with high confidence 
levels may be increased significantly. 

The only advantage that TM imagery appears to have over MSS 
data within this study is the ability to detect smaller canopies 
and more clearly delineate larger canopies. The higher spatial 
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resolution of TM Imagery could further Increase In Importance for 
study areas located In Upper Chesapeake Bay where the SAV canopies 
tend to be smaller and less numerous than In the Lower Chesapeake 
Bay. 


CONCLUSIONS 

Landsat TM and MSS imagery are similarly effective in de- 
tecting SAV in Lower Chesapeake Bay occupying bright, sandy sub- 
strates. The increased radiometric sensitivity and spectral reso- 
lution and the addition of a third optical band does not afford TM 
imagery a significant advantage over MSS data. However, the reader 
is advised against extending this conclusion to other topics in 
water research or, for that matter, SAV occupying substrates having 
a spectral reflectance more closely resembling the canopy. For 
cases in which spectrally similar water masses are to be delineated 
or when SAV is to be detected upon dark, muddy substrates, the 
radiometric and spectral characteristics of the TM bands may offer 
a significant improvement over those of the MSS. 

The only apparent advantage to using TM imagery over MSS 
data for detecting SAV in Lower Chesapeake Bay is that, by virtue 
of higher spatial resolution, smaller SAV canopies may be detected. 


in 
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Table 1. 

Radiometric Sensitivity* Of Optical TM And MSS Bands 
Wavelength Region 


Band 

(nm) 

Gain 

Offset 

MSSl 

500-600 

5.570 

-1 .1 1 u 

MSS 2 

600-700 

7.21 6 

-2.886 

TMi 

i<50-520 

15.9*111 

2. *123 

TM2 

530-610 

8.199 

2.329 

TM3 

620-690 

10.81*1 

1.265 


1 COUNT - 

GAIN * L ♦ OFFSET 



COUNT - 

0-255 



GAIN - 

OOUNT / mwcnT * y m" * s tr “ 

1 


L * 

SCENE RADIANCE 



OFFSET - 

0-255 



Table 2. 

Maximum Cluster Radius (R) And Minimum Between-Cluster 
Distance (D) Used To Classi. y The TM and MSS Imagery 



TM 

MSS 

R 

.903 

.580 

D 

1.806 

1 .1 60 
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ORIGINAL PAGE IS 

Of POOR QUALITY 


1 SAV 
Crown 
Cover 

0-10 

10-H0 

U0-70 

70-100 


1 SAV 
Crown 
Cover 

0-10 
10-110 
HO- 70 
70-100 


Table 3. 

Landaat TM and MSS Classification 
Accuracies For All Water Areas 


Confidence 


Sensor 

0-251 

25-501 

50-751 

— J 
i 

o 

o 

TM 

100.0 

0.0 

0.0 

0.0 

MSS 

97.3 

2.7 

0.0 

0.0 

TM 

72.1 

17.5 

9.2 

1.2 

MSS 

76.2 

8.3 

15.5 

0.0 

TM 

100.0 

0.0 

0.0 

0.0 

MSS 

92. n 

7.6 

0.0 

0.0 

TM 

98.5 

1.2 

0.1 

0.2 

MSS 

96.7 

1.8 

1.0 

0.5 


Table H. 

Landaat TM and MSS Classification Accuracies 
For Water Areas Less Than 6 Feet 


Confidence 


Sensor 

0-251 

25-501 

50-751 

75-1001 

TM 

100.0 

0.0 

0.0 

0.0 

MSS 

97.3 

2.7 

0.0 

0.0 

TM 

72.1 

17.5 

9.2 

1.2 

MSS 

76.2 

8.3 

15.5 

0.0 

TM 

99.9 

0. 1 

0.0 

0.0 

MSS 

82.9 

17.1 

0.0 

0.0 

TM 

53.6 

16.7 

10.0 

19.7 

MSS 

55.9 

15.7 

15.1 

13.3 
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1. INTRODUCTION 

This research Is concerned with assessing the effectiveness of 
LA.1DSAT TM and MSS sensors for detecting submerged aquatic vegetation 
(SAV). Me approach the problem from a theoretical standpoint In which 
we simulate volume reflectance from SAV under a variety of conditions. 

Here, we focus upon the spectral and radiometric qualities of TM bands 

1, 2 and 3 and MSS bands 4 and S. 

2. BACKGROUND 

The vast majority of remote sensing research Involved with detecting 
surface features addresses the problem In a purely stochastic fashion. 
Typically, computers are used to generate categories within digital 
multi-band data based upon spectral signatures. Signatures may be 
derived either manually by selecting subsets of data that correlate with 
ground observations or by Instructing the computer to group similar 
signatures. The Image Is then classified by statistically assigning 
each pixel to one of the categories. These techniques are nonparametrlc 
In the sense that the classifier makes no assumptions about the physical 
characteristics of surf act features. Under certain conditions, e.g. 
large, homogeneous surface features relative to the spatial resolution 
of the sensor, the nonparametrlc approach has been very successful . 

Numerous examples are available In the literature and Molk (1980) gives 
a good overview of the subject. 

In a parametric approach, the researcher attempts to predict upwelllng 
radiance at an airborne or orbiting sensor by modeling radiative transfer. 
The problem becomes significantly more complicated because the researcher 
requires knowledge of the t«?nsm1ttance and reflectance characteristics 
of the surface feature as well as the atmosphere. The advantage of the 
parametric approach Is that the researcher gains valuable insight for 
Interpreting observed radiance changes in the Imagery. 

3. METHODOLOGY 

3.1 Spectral Assessment 

Radiative transfer theory Is used to model upwelllng radiance that 
would be received by LANDSAT TM and ICS viewing a Hypothetical estuarine 
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environment, shown In Figure 1. The environment Is composed of e deer 
maritime atmosphere, an optically shallow estuary of either dear or 
turbid water, and three possible bottom types: vegetation, sand or mud. 

3.1.1 The Atmosphere 

A clear maritime air mass Is selected and Is described by fiuttman 
(1968). Using solar Irradlance data (Bast et el., 1965) es Input, solar 
and sky Irradlance Is calculated at sea level as well as path radiance 
that would be received by an orbiting sensor. These calculations assume 
a solar zenith angle of 50°. 

3.1.2 The Water Column 

Irradlance reflectance of the water column and bottom Is calculated 
using a quasi-single scattering model developed by Phllpot (1981). 

v£o-#-V)*v* v <’> 

* . 

Where Rj • Irradlance reflectance of the water and bottom, 

« Irradlance backscatter coefficient, 

IL ■ total Irradlance attenuation coefficient (the sum of the 

T attenuation coefficients for upwelllng and downwelllng 
irradlance), 

d - water depth, 

A^ ■ Irradlance reflectance of the bottom. 

Equation (1) Is used, at the expense of slightly lower accuracy, 
rather than a two-flow or Monie Carlo approach because the two necessary 
terms, B, and tL, are easily calculated from simple field measurements 
available In the literature and the amount of computer time required 
Is relatively small. 

The water column Is assumed to be vertically homogeneous and the 
water surface calm and flat. Two very different water qualities are 
considered; clear ocean water and turbid fresh water. Measurements 
representing clear oceanic water are documented In Tyler et al. (1972). 
Measurements representing very turbid water were made In San Vincente, 
a man-made lake northeast of San Diego, and are documented in Tyler 
and Smith (1970). Calculated values of 8. and K t for both water types 
are shown in Table 1. 

3.1.3 The Bottom Types 

Three different reflectance profiles are selected from the 
literature to represent a submerged plant canopy, sand and mud. The 
spectral reflectance of each bottom type Is shown In Figure 2. 
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3.1.4 Upwel ling Radiance Above the Earth 

The upwel ling radiance that would be received by an orbiting sensor 
Is calculated as: 

L o $ -l H o R l ♦ V l p • (2 > 

Lq$ » upwel ling radiance at the orbiting sensor, 

H * the combined effects of the atmosphere and air/water interface 
0 upon radiance reflected from the water, 

Rj » Irradiance reflectance of the water and bottom, 

P ■ Fresnel reflectance of downwelllng radiance at the water surface, 

fl 

L^- « downwelllng sky radiance, 

T 0 • optical thickness of the atmosphere in the zenith direction, and 
Lp « atmospheric path radiance received by the sensor. 

Upwel ling radiance values are calculated every 10 nm over the wave- 
length range 450 nm to 700 nm. Radiance values are then summed over 
TM bands 1 (450-520 nm), 2 (520-600 nm), and 3 (630-690 nm) and MSS bands 
4 (500-600 nm) and 5 (600-700 nm). 

3.1.5 Spectral Quality Index 

A spectral quality Index Is defined based upon the work of Lyzenga 
and Polcyn (1978) as: 

SQ1 « |L t - L b I L b ‘ 1/2 (3) 


where: 

Lj • upwel ling radiance representing the submerged plant canopy, and 

L b * upwel ling radiance representing a background of sand or mud. 

Equation (3) is quite similar to the equation for apparent contrast. 
Relative values of SQI between bands indicate relative effectiveness 
In detecting submerged plant canopies. 

3.2 Radiometric Assessment 

The radiometric resolutions of TM bands 1, 2 and 3 and MSS bands 
4 and 5 are compared for the same cases discussed In the previous section. 
Upwel ling radiance at satellite altitude calculated for each band Is 
converted to count values using calibrations taken from Barker et el. 
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(1983) for TM sensors and Bauer (1980) for LANDSAT-3 MSS sensors 
(see Table 2 for conversion constants). Apparent contrast Is then 
calculated between target and background signals as: 

AC * |Cy * Cg| 

where 

Cy • count value calculated from target signal, and 

Cg * count value calculated from background signal. 

3.3 Quick-Look Evaluation of TM Data 

An EROAS 400 microcomputer Is used to enhance and classify a TM 
Image of Broad Creek, Maryland. The Image was generated on November 2, 
1982. The scene Is classified using an unsupervised clustering routine. 
The classified scene Is compared with aerial photography of the area 
collected on August 2, 1982. Of primary Interest In the scene Is a 
medium size SAV bed located on the west bank of Broad Creek, latitude 
38°44'20 M N, longitude 76°15'25"W. The SAV was photographed from a 
Cessna 174 at low altitude using a hand-held camera and Kodak Pan-X 
black and white film.. 

4. RESULTS - 

SQI values are calculated for clear water overlying a sand bottom, 
clear water overlying a mud bottom, turbid water overlying a sand bottom, 
and turbid water overlying a mud bottom. The results are shown In 
Tables 3-5 respectively. In all cases, when the water depth Is optically 
shallow TM band 3 and MSS band 5 are most effective In detecting the 
submerged vegetation. This Is Intuitively correct since the Inherent 
contrast between the vegetation and both the sand and mud Is a maximum 
In these two bands. In clear water, as the water depth Increases, the 
optimum band shifts from TM band 3 and MSS band 5 to TM band 1 for 
both bottom types. Here, the attenuation of bottom signal is less 
In TM band 1 than in TM band 3 and MSS band 5. In turbid water the 
relative effectiveness of each band is similar for all depths considered. 

Spectrally, MSS band 4 Is quite similar to TM band 2 and MSS band 5 
Is similar to TM band 3. Nevertheless, some unanticipated changes In 
relative SQI values between these bands are observed with respect to 
Increasing water depth. Shown In Table 5, MSS band 4 Is slightly more 
effective than TM band 2 in discriminating between submerged vegetation 
and sand at 0.5 meters. When the water depth increases to 1.5 meters, 

TM band 2 becomes slightly more effective. Similar results are seen 
In Tables 2 and 3 between TM bands 3 and MSS band 5. 

In most cases, SQI values decrease with Increasing depth. This 
also seems Intuitively correct as the apparent contrast decreases with 
either Increasing attentuatlon or Increasing pathlength. In Table 4, 

SQI decreases In TM band 2 between 0.5 and 1 . : meters and then Increases 
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between 1.5 meters and 10.0 meters. To understand what Is happening here 
It Is necessary to know something about the variation In the optical charac- 
teristics of the water and the bottom types across the band. Figure 3 
shows the spectral reflectance of the vegetation and mud within TM band 2. 
Whereas the largest reflectance from the vegetation occurs In the shorter 
wavelength portion of the band, the reflectance of the mud Increases toward 
the longer wavelength region. The average reflectance of the mud Is 
slightly larger than the vegetation so that In the absence of a water 
column, the mud would appear brighter. Also shown In Figure 3 is the total 
attenuation coefficient of the clear ocean water which increases signifi- 
cantly toward the longer wavelength region of the band. Under these 
circumstances, the signal from the mud decreases rapidly with Increasing 
water depth than does the vegetation signal. Figure 4 is a plot of SQI 
without taking the absolute value of the numerator. Values were calculated 
for several depths between 0 and 20 meters. In very shallow water, the 
signal from the mud Is greater than the signal from the vegetation and the 
Index takes on a negative value. At some Intermediate depth slightly greater 
than 1.0 meters, the two signals are equal and SQI Is zero. At still 
greater depths the vegetation appears brighter than the mud and SQI values 
are positive. As the depth becomes very deep, both signals take on the 
vilue of optically deep water and SQI again falls to zero. This emphasizes 
the Importance of making optical measurements within natural waters at 
appropriate spectral resolutions. If all optical measurements were made 
with a bread band radiometer representing TM band 2, the intermediate zero 
contrast would never have been noticed and the model predictions of up- 
welling radiance would have been In gross error. 

Also shown in Tables 3-6 are apparent contrast calculations In which 
upwelllng radiance Is converted to count values. The results simulate the 
absolute count difference within the satellite Image between pixels repre- 
senting the SAV and those representing the surrounding unvegetated sand or 
mud. For the case of clear ocean water, TM band 1 offers the best discri- 
mination between the SAV and the background, both sand and mud. This 
result does not agree with that of the SQI analysis, the reason being that 
each sensor has a unique counts-to-radiance calibration. In Table 2, the 
higher the sensor gain, the more sensitive It Is to changes In reflectance. 
Thus, even though spectrally TM band 3 Is better than TM band 1, the larger 
gain of TM band 1 more than makes up for the difference. This same effect 
occurs between TM bands 2 and 3 and the spectrally similar MSS bands 4 and 
5. In each case the gain Is larger for the TM bands, making them signifi- 
cantly more efficient In detecting the SAV. 

A TM band 5 image of Broad Creek, Maryland is shown in Figure 5a. 

The image was generated on November 2, 1933. Figure 5b shows an aerial 
photograph of a SAV bed located along the west shore of Broad Creek. The 
SAV appears as a dark area adjacent to the shoreline on an otht'"wise bright 
bottom. The photo was taken on August 2, 1982, three months pno”* to the 
TM Image. The TM Image is classified with an unsupervised clustering 
routine using an ERDAS 400 microcomputer. The classified Image Is shown 
In Figure 6. The area classified as SAV correlates well with that shown 
in the aerial photography. . 
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5. CONCLUSIONS 


The effectiveness of an orbiting sensor In discriminating spectrally 
between submerged features Is a function of the Inherent contrast between 
the submerged features and how strongly the bottom signal Is attenuated 
by the water column. In optically shallow water the Inherent contrast 
Is the controlling factor. Thus* the optimum sensor band Is that which 
correlates with the greatest Inherent contrast between the submerged 
features. In optically deeper water* the optimum sensor band Is that 
for which the water column Is most transparent. 

In the clear ocean water the optimum band for detecting vegetation 
on a sand or mud bottom Is shown to change with the optical depth of the 
water. In the turbid San Vicente water the optimum band In optically 
shallow water remains the optimum band In optically deeper water. 

Under certain conditions the spectral quality Index will decrease 
to zero for some intermediate depth and then Increase for yet deeper 
depths. This could have a significant Impact upon how changes in water 
reflectance should be Interpreted and can only be predicted with detailed 
knowledge of the spectral variation 1 of optical parameters across the 
sensor band. 

Radiometric calibration is shown to be a dominant factor In deter- 
mining which TM or MSS band will be optimum for detecting a submerged 
feature such as SAV. Sensors with higher gain settings In general 
produce data with more Information. The gain built Into the TM sensors 
are significantly higher than those of the MSS, giving TM data a -clear 
advantage for detecting submerged features. In several cases, a band 
Is found to be optimum for detecting SAV In a radiometric sense and yet 
sub-optimal In a spectral context. 

The Thematic Mapper appears to be useful for detecting SAV in the 
case of Broad Creek, Maryland. A nonparametric classification of a 
November 2, 1982 image clearly discriminated between SAV and the sur- 
rounding unvegetated bottom. 
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TABLE 2 


Counts/Radiance Conversion Constants for LANDSAT-4 TM bands 1, 2 and 3 
(Barker et al., 1983) and LANDSAT-3 MSS bands 4 and 5 (Bauer, 1980). 


Band 

Dynamic Range 

R m1n* ** 



(mw/cm -str-um) 

^nax 

Gain 

2 

(counts/mw/(cm - 

TM 1 

-0.152 

15.842 

15.78 

TM 2 

-0.284 

30.817 

8.10 

TM 3 

-0.117 

23.463 

10.62 

MSS 4 

0.4 

22.0 

5.79 

MSS 5 

0.3 

17.5 

7.27 


* Largest radiance value observed for a count value of 0. 

** Smallest radiance value observed for a TM count value at 255 or an MSS 
count value of 126. 
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Table 3 


SQI VALUES* AND APPARENT CONTRAST** BETWEEN SAV AND SAND FOR 
THE CASE OF CLEAR OCEAN WATER 


Band 

d s 0.5m 

d*l. 

5m 

d»10.0m 

351 

A C 

351 

aC 

351 

IT 

TM 1 

143.8 

16.3 

128.9 

14.4 

49.2 

5.2 

TM 2 

101.4 

5.6 

79.9 

4.3 • 

12.5 

0.6 

TM 3 

235.8 

14.0 

120.7 

6.0 

0.1 

0.0 

MSS 4 

108.2 

3.7 

88.1 

2.9 

18.6 

0.6 

MSS 5 

214.1 

8.2 

117.4 

3.8 

0.5 

0.0 


* Units ■ (mw/cn^-str-nm) 1 ^ x 10'^ 
** Units » Counts 


Table 4 

SQI VALUES* AND APPARENT CONTRAST** BETWEEN SAV AND MUD FOR 
THE CASE OF CLEAR OCEAN WATER 


d*0.5m d=1 . 5m d*10.Qm 


Band 

SQI 

AC 

SQI 

AC 

SQI 

AC 

TM 1 

70.5 

7.5 

62.9 

6.7 

23.4 

2.4 

TM 2 

2.8 

0.4 

1.8 

0.1 

4.2 

0.1 

TM 3 

151.8 

7.9 

74.5 

3.4 

0.01 

0.0 

MSS 4 

13.23 

.6 

8.5 

0.5 

0.2 

0.1 

MSS 5 

127.5 

4.3 

67.9 

2.0 

0.2 

0.0 

* Units 

2 

« (mw/cm - 

str-nm)^ 

x 10° 





** Units ■ Counts 
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Tabic S 


SQI VALUES* AND APPARENT CONTRAST** BETWEEN SAV AND SAND 
FOR THE CASE OF TURBIO SAN VICENTE WATER 


Band 

d»0.5m 

d-1 

• 5m 

d* 

10.0m 


5T 

SOI 

&c 


AC 

TM 1 

69.4 

7.3 

31.5 

1.4 

14.4 

0.0 

TM 2 

70: 6 

3.7 

42.8 

1.2 

25.6 

0.0 

TM 3 

174.3 

9.9 

88.9 

2.1 

43.3 

o.c 

MSS 4 

71.7 

2.3 

42.1 

0.7 

24.5 

0.0 

MSS 5 

154.8 

5.8 

83.1 

1.3 

42.1 

0.0 

* Units 

2 

» (mw/cm • 

-str-nm)^ 

x 10‘ 3 




** Units 

» Counts 







Table 6 * 

SQI VALUES* AND APPARENT CONTRAST** BETWEEN SAV AND MUD 
FOR THE CASE OF TURBID SAN VICENTE WATER 


Band 

d»0 

.5m 

d*l 

.5m 

d* 

10.0m 

SOI 

AC 

SQI 

AC 

SQI 

AC 

TM 1 

70.5 

3.3 

62.9 

0.6 

23.4 

0.0 

TM 2 

3.1 

0.3 

1.6 

0.1 

0.8 

0.0 

TM 3 

109.0 

5.6 

53.5 

1.2 

25.5 

0.0 

MSS 4 

8.5 

0.4 

4.3 

0.1 

2.2 

0.0 

MSS 5 

91.4 

3.0 

46.7 

0.7 

23.3 

0.0 

* Units 

« (mw/cm - 

str-nm) 1 ^ 

x 10‘ 3 





** Units ■ Counts 









450 500 550 600 650 700 

Wavelength (nm) 


Figure 2. Xrradiance reflectance of three bottoa types used in 
the water reflectance model (Philpot, 1981) . Sand ■ 
beach sand (Lyzenga and Thomson, 1978) , mud ■ dark 
soil (Lyzenga and Thomson, 1978) , and vegetation * an 
infinitely deep Ilex canopy (Gates et ml, 1965). 
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Figure 4. SQI values calculated without the absolute value of 

the nuserator for TM band 2. The water type is clear 
ocean, L is the signal representing subaerged vege- 
tation, and Lg is the background signal representing 
a bottoa of aud. 
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FIGURE 5. 

A portion of Chesapeake Bay it shown in plate i and an enlarged 
portion of Bruec Creek. Maryland appears in plate b. The boxed 
area within plate b was photographed from low altitude on 
Augurt 2. 1982. The photograph, plate c, shows a medium 
size SAV bed adjacent to the shoreline. 
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FIGURE 6. 


A port] or of the November n , 1382 TM image clashed using an ERDAS 400 micro- 
computer and an unsupervtsed clustering routine. The cashed lines represent the area 
photographed from low altitude and shown in Figure 5c. Pixels classified as SAV are 
shown as I. optically deep water as Q submerged sand as □, and land as either □ 
(fields) or E3 (trees). 
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ABSTRACT 

Spectral data, simulating thematic mapper bands 3, 4 and 5 were 
gathered In salt and brackHh marshes using a hand-held radiometer. 

Simple regression models were developed equating spectral radiance 
Indicias with total Idve biomass for S. alternl flora In a salt marsh ^nd 
for a variety of plant species in a brackish marsh. Models were then 
tested using an Independent set of data and compared to harvest estimates 
of biomass. In the salt marsh, biomass estimates from spectral data 
were similar to harvest biomass estimates during moSt of the growing 
season. Estimates of annual net aerial primary productivity calculated 
from spectral data were within 211 of production estimated from harvest 
data. During August, biomass estimates from spectral data in the brackish 
marsh were similar to biomass estimated by harvesting techniques. At 
other times during the growing season, spectral data estimates of biomass 
were not always comparable to harvest biomass estimates. Reasonable 
estimates of wetlands biomass are possible during the peak of the growing 
season (August) using spectral data similar to thematic mapper bands 3, 

4 and 5 gathered with hand-held radiometers. 

INTRODUCTION 

Basic to the understanding of wetland function and value is the 
quantification of energy fixation. Reduced carbon compounds comprising 
macrophytic biomass provide the energy necessary to maintain the hetero- 
trophlc organisms which feed upon them. Aboveground biomass represents 
only a portion of total net primary production (belowground production 
can also be substantial) however, the fixed carbon In this biomass Is 
the surplus available to heterotrophlc organisms and Is readily detected 
with remote sensing devices. 


Presented at the Landsat-4 Science Characterization Fany Results 
Symposium, NASA Goddard Space Flight Center, Greenbelt, MD, 23-24 
February 1984. NASA Conference Publ . I V : 251-269 . 


Many salt marshes on the eastern coast of the J.S. are dominated by 
a single plant Spartlna alternlflora , Lolsel. (Relmold 1977). _S. 
alternlflora biomass can vary from near 100 to over 3000 gdw nf^dependlng 
upon substrate type, soil salinity. Inundation frequency or other edaphlc 
factors. As Interstitial water salinity decreases to the 18-15 o / 0O 
range, S. alterni flo ra Is gradually replaced by a variety of plant 
species. In genera), the lower the salinity, the greater the diversity 
of the resident plant population. The monospecific gramineous canopies 
of the salt marsh are transformed to brackish marsh canopies comprised 
of mixtures of gramineous, broadleaf and leafless plants. Spatially and 
temporally the canopies encountered In the brackish system are complex 
and dynamic. The morphologic diversity of the plant conmuni ty produces 
an equally diverse spectral signature. 

The study described herein was designed to develop simple regression 
models equating spectral radiance Indlcles with plant biomass. The 
radiance data were spectrally similar to thematic mapper bands 3, 4 and 
5 and have been suggested as being superior to MSS wavebands for vegetation 
monitoring (Tucker 1978). Our first objective was to use these models 
to predict biomass and net primary productivity for a sale marsh and to 
predict biomass for a brackish marsh using '” r ound gath;red spectral 
radiance data. Our second objective was o ess the validity of our 
biomass estimation by comparing the spectrally estimated biomass values 
with biomass estimates obtained by traditional harvest techniques. 

METHODS 

Salt Marsh 

A portion of the Canary Creek salt marsh in Lewes Delaware was 
selected for study (Figure 1). The marsh was dominated by monospecific 
stands of S. alterniflora with relatively small areas of Pi s t i chi i s 
spicata (LTj Greene . Occassional ly, Sal i corn i a europaea I~6r Limonium 
sp. occurred mixed with S. alterni florae Six transects extending from 
the creek edge to the upland were established. The transects were 
spaced approximately 210 meters apart and stations were designated every 
30m along each transect. A total of 40 stations for the whole marsh 
were established. By establishing stations in a systematic manner along 
transects, we sampled a representative cross section of S. alternlflora 
height forms occurring within the marsh. Beginning on l5~May 198T and ~ 
continuing every three weeks until 9 OctoDer 1981, four stations were 
selected, using a table of random numbers, from each transect. For each 
of the eight sampling dates, 24 stations were sampled yielding a total 
of 192 for the season. 

Brackish Marsh 


A brackish marsh near the headwaters of Old Mill Creek at Lewes, 
Delaware was selected for study (Figure 1). The Interstitial water 
salinity varied from approximately 15-18 u / 00 at the downstream end to 
about 12-10°/oo at the upstream end. Plant species comprising the plant 
comnunlty changed with apparent soil salinity. Four transects extending 
from the creek edge to the upland were established 180 meters apart 
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along the long axis of the creek. Stations were spaced at 30m Intervals 
along each transect. A total of 29 stations were designated In the 
marsh. The marsh was sampled 3 times during the growing season (June, 
August, September) at approximately 6 week Intervals. At each sampling 
date, 4 stations were selected using a table of random numbers from each 
transect yielding 16 stations and a season total of 48. 

Spectral Radiance Data Collection 

Before harvesting, a GSFC Mark II hand-held radiometer (Tucker et 
al. 1981a) was used to measure canopy radiance directly over each area 
to be harvested. The radiometer contained three wavebands spectrally 
configured with interchangeable interference filters to match bands 3, 4 
and 5 of the Landsat-4 thematic mapper. The sensor head housed a red 
band (0.63 - 0.69 ym, RED) sensitive to chlorophyll concentration, a 
near Infrared band (0.76 - 0.90 ym, NIR) sensitive to plant tissue 
structure or biomass and a middle infrared band (1.55 - 1.75 ym, IR) 
sensitive to leaf moisture. Data were recorded simultaneously for all 
three bands. 

No more than 5 days prior to the actual harvest, spectral radiance 
was determined for each plot. The radiometer was leveled approximately 
1.5 meters above the top of the plant canopy. In the brackish marsh a 
wooden step ladder was employed to achieve proper Instrument height. 

The radiance was measured 3 times over each plot. Radiance data were 
not collected under cloudy or very windy conditions and were always 
collected within 2 hours of solar noon. Radiance data were collected 
during -low tide, however, there were occasions when small amounts of 
tidal water remained pooled on the marsh. Standing water was avoided 
whenever possible due to potential specular reflectance from the water 
surface. 


Spectral radiance data were transformed and expressed as a normalized 
difference of two bands as outlined by Krlegler et al. (1969) and Rouse 
et al . (1973). The red and near Infrared radiance values were combined 
in the following manner: 


VI 


NIR - RED 

jrnr+m 


where VI Is the vegetation Index, NIR is the near Infrared band radiance 
and RED Is the red band radiance. A similar combination of the near 
infrared and middle Infrared bands was performed (middle infrared substituted 
for RED in the above expression) and termed the infrared Index (II). 

The infrared index is used here with the understanding that the near 
infrared and the middle infrared bands are spectrally different and that 
normalization in this manner may not be totally valid since adjacent, 
spectrally similar bands were assumed for this transformation (Kriegler 
et al . 1969). Index values were preferred to raw radiance data because 
the normalization procedure tends to compensate for changes in solar 
irradlance caused by seasonal changes In solar zenith angle and/or 
atmospheric conditions (Tucker et al . 1979a). 
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Harvesting Procedures 

After collection of radiance data, all vegetation (Including standing 
dead material) within a 0.25 m 2 frame was clipped at soil level at each 
station. The plant material was bagged and returned to the laboratory 
for processing. A 1/3 subsample (by wet weight) was drawn from the salt 
marsh samples and sorted Into live and dead components. The brackish 
marsh samples were sorted completely by species, then only plant species 
with a large volume of material were subsampled and finally live and 
dead components determined for each species. Host plants In the brackish 
marsh samples contained little attached dead so the majority of the dead 
material was considered as a composite of the whole sample. The subsampling 
procedure for both marsh types greatly reduced processing time. We felt 
that a representative subsample was sufficient to determine relative 
proportions of live and dead tissue and by retaining the entire 0.25 m 2 
sample for biomass determination, we preserved the best estimate of 
biomass, given the oftentimes spatially heterogeneous distribution of 
biomass within the area sampled. Proportions of live and dead tissue 
determined from the subsample were then extrapolated to the entire 
biomass sample. All plant material was dried at 60°C to a constant 
weight, weighed to the nearest O.lg and expressed as grams dry weight 
per square meter (gdw m’ 2 ). 

RESULTS 

Salt Harsh - Biomass Estimation 

Regression models equating S. alterni fl ora biomass and spectral 
radiance indicies were developed~during the 1980 growing season (Hardisky 
et al. 1983a). These models Included short and tall form $. alternjflora 
sampled from June through November using the hand-held radiometer. 

Table 1 lists the models for the vegetation and Infrared indicies. To 
estimate live and dead biomass, one solves for live .leaf biomass (L) in 
the first equation, for total live biomass (live leaves and stems, T) in 
the second equation and then substitutes the L value (from equation 1) 
into the last equation solving for dead biomass (D). Through this 
sequence, the necessary parameters for annual net aerial primary produc- 
tivity estimation (total live biomass and dead biomass) are calculated. 

Figure 2 depicts biomass estimates from harvesting and computed 
from the vegetation index for S. alterni flora . For live and total 
biomass, the predicted values were similar to the harvest values throughout 
the year. Live biomass estimates using the vegetation index ter J 2 d to 
be lower than harvest estimates during the first half of the growing 
season and were higher than harvest estimates during the latter half of 
the growing season. Peak biomass was attained in early August at which 
time canopy development would be the fullest. We would expect the 
greatest proportion of live biomass conpared to dead biomass to occur at 
this time and for the greatest amount of live tissue to be in the upper 
portions of the canopy. This may have contributed to the higher vegetation 
Index biomass estimates seen at peak biomass and thereafter. Dead 
biomass did not vary greatly over the season (from 260 to 400 gdw m" 2 ) 
however, decomposition of carry-over dead material (from previous growing 
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season, high In the canopy) and the concurrent addition of dead material 
from Immature culm mortality and leaf senescence (low In the canopy) can 
change the relative vertical position of the dead biomass without large 
changes In biomass thereby altering reflectance from this component. It 
would seen reasonable that this change In the location of dead material 
contributed to the lower live biomass predictions In the first part of 
the season (when more dead material was In the canopy) and to the higher 
live biomass predictions in the latter part of the season (when less 
dead material was In the canopy). 

Table 1 

REGRESSION MODELS FOR PREDICTING S. ALTERN I FLORA BIOMASS 

Coefficient 
of Determination 


Regression M odel (r 2 ) 

VI = .382 + .068 1 n ( L ) 0.75 
VI = .149 + .096 ln(T) 0.64 
VI = .760 + .055 ln(L/D) 0.88 
II * .178 + .104 ln(L). 0.87 
II =-.185 + .148 ln(T) 0.76 
II = .752 - .078 1 n( L/D) 0.88 


L - Live leaf biomass: T - Total live biomass; D - Dead biomass 
Units on all biomass = grams dry weight per meter squared 
n for all regression models * 96 
VI = Vegetation Index; II = Infrared Index 

Harvest and infrared index estimates of live and total S^ al terni flora 
biomass are compared in Figure 3. The infrared index estimates of live 
biomass were very close to the harvest estimates during the early part 
of the growing season. After peak biomass, the infrared index estimates 
of live biomass were lower than harvest estimates. The reasons for the 
apparent underestimates at the end of the growing season are not clear, 
however, we would postulate that the lower water content of the plant 
tissue due to tissue maturity and increased interstitial water salinity 
at this time in the growing season would cause an increase in the middle 
infrared reflectance resulting in lower infrared index values. Overall 
the infrared index and vegetation index estimates of biomass were similar. 

Although the mean biomass estimates from harvesting and from radiance 
indicles were similar, considerable variability existed for Individual 
measurements during some parts of the growing season. Figure 4 presents 
coefficients of determination (r 2 ) between harvest estimates and vegetation 
or Infrared Index estimates of live (a) and total (b) biomass. High r 2 
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values indicate good linear association between Individual harvest and 
radiance Index estimates of biomass. Low r 2 values Indicate considerable 
variation existed for some samples between the two biomass estimates. 

From May through June, the most disagreement between harvest and radiance 
Index biomass estimates existed. August was the best time period for 
agreement of the biomass estimates with the latter part of the growing 
season as a whole being better than the first part. The harvest and 
radiance Index estimates of live biomass were always In better agreement 
than the total biomass estimates. Harvest and radiance Index biomass 
means In Figures 2 and 3 were generally very similar. This suggests 
that during those periods when low r 2 values existed between harvested 
and predicted biomass, the amount of overestimation and underestimation 
were similar among the 24 samples, yielding means which were comparable. 
When using radiance index values for biomass prediction, one should 
scrutinize the results from Individual measurements to assure that each 
estimate Is reasonable, particularly early In the growing season. 



Fig. 2 Vegetation Index and harvest estimates of S. alterniflora 
biomass over the growing season. Bars representing 1 standard 
error of the mean are shown whenever they exceed the size of the 
point symbol. Points are slightly offset to Improve readability. 
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Fig. 3 Infrared Index and harvest estimates of S^ alterniflora 
biomass over the growing season. Bars representing 1 standard 
error of the mean are shown whenever they exceed the size of the 
point symbol. Points are slightly offset to improve readability. 

Salt Marsh - Productivity Estimation 

Using the mean biomass estimates shown In Figures 2 and 3, annual 
net primary productivity was computed for the S. al terni flora salt 
marsh. We employed established production calculation techniques used 
widely in salt marsh systems (Table 2). There is some discussion as to 
the adequacy of these techniques and we recognize the shortcomings of 
each. However, these techniques are widely used anu will yield values 
comparable to other studies. Linthurst and Reimold (1978) provide an 
excellent comparison of various harvest techniques for estimating net 
aerial primary productivity in estuarine marsh systems. 

Annual net aerial primary productivity estimates for the Canary 
Creek marsh are found in Table 2. The vegetation Index estimate of 
productivity was within 4? of the Smalley (1958) harvest technique and 
no more than 21i different fror the other techniques. The Infrared 
Index estimates of productivity were within 205 of the Smalley (1958) 
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technique and were 4-5* different from the other techniques. The key to 
the agreement between the productivity estimates computed from harvest 
and radiance Index data was the high correlation between predicted and 
harvest biomass during peak biomass (August). For all but the peak 
standing crop method of production estimation, successive positive 
Increments In biomass are sumned throughout the growing season, normally 
culminating at the point of peak biomass. Therefore, the estimation of 
peak biomass Is the single most critical measurement In determining 
annual production. 


Table 2 

annual net aerial primary productivity estimates 

FOR SPARTINA ALT ERNIFLORA IN CANARY CREEK MARSH 

Remote Sensing Estimate 

Method Harvest Estimate VI II 


peak standing crop 

517 

600 

489 

Milner & Hughes (1968) 

523 

661 

498 

Morgan (1961) 

517 

600 

497 

Smalley (1958) 

634 

661 

506 


all values are grams dry weight per square meter per year 
VI * Vegetation Index; II * Infrared Index 

Brackish Marsh - Model Development 

Modeling the relationship between spectral radiance indlcles and 
live aerial biomass in brackish marsh plant coumunltles required consid- 
eration of diverse morphologic characteristics among plants residing in 
the same cormunity. To Illustrate this point. Figure 5 shows the linear 
relationship between vegetation Index and live biomass for a variety of 
wetland plants. Iva , Polygonum , and Solldago represent broadleaf or 
deciduous canopies and exhibit rapid Increases In vegetation Index for 
relatively small changes In biomass. This characteristic suggests that 
the spectral Index can become saturated rapidly. Structurally, these 
canopies maintain most leaf surfaces In the horizontal plane and generally 
form a complete canopy cover which reduces or eliminates exposure of 
dead components or soil background to solar irradiance. This combination 
of canopy characteristics yields a very absorptive canopy in the red 
region and a very reflective canopy In the near infrared region, thus 
the high vegetation index relative to the amount of live biomass present. 

The opposite extreme to the broadleaf canopy would be the leafless 
canopies represented by Sal Icornla and Sclrpus . Both Sal Icornla vlrglnlca 
and Sclrpus olneyl possess erect, ^leafless stems with most green tissue 
In the vertical pTane and primarily soil background and dead plant material 
In the horizontal plane. Normally these canopies are very open with soil 
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surface characteristics potentially contributing greatly to the 
observed spectral radiance. 

The third canopy type represented In Figure 5 Is the gramineous canopy 
type of Spartlna and Typha . Spartlna alternlflora exhibits alternate 
leaves along the 1 ength of the stem whereas Typha a noustlfolla has basal 
leaves. Both plants form canopies with portions o T leaves in the horizontal 
and In the vertical plane. S^ alternlflora exhibits a broad range of 
canopy configurations as a result of Its wide environmental tolerance 
limits. Both plant canopies can maintain substantial quantities of dead 
material within the canopy. The amount of live leaf tissue determines 
to what degree dead material and soil background will Influence spectral 
reflectance. Theoretically, the occurrence of flat leaves (portions of 
which may be horizontal) In the gramineous canopy would place them 
somewhere between the broadleaf and leafless canopies In terms of an 
Increase In vegetation Index value for an Increase In biomass (l.e. an 
Intermediate slope). In practice this does not occur because the measured 
vegetation Index represents the composite of reflectance from vegetation 
(live and dead) and the soil. In the case of the gramineous canopy, the 
dead vegetation and soil are oftentimes well Illuminated and contribute 
significantly to the measured vegetation Index. The net effect Is a 
lessening of vegetation Index Increases with Increasing biomass. 



Fig. 5 Comparison of the relationship between vegetation Index values 
and amounts of live biomass for a variety of wetland plants. Iva, 
Polygonum and Sol Idago represent broadleaf canopies, Typha and Spartln a 
(dashed lines) represent gramineous canopies and Sallcornla and 
Scirpus represent leafless canopies. 
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Brackish marsh plant convnunltles are often composed of numerous 
plant species with morphologies represented by each of the groups discussed. 
The problem of equating spectral data with live biomass becomes a function 
of the species composition of the community In question. At present, we 
lack sufficient data for modeling specific mixtures of the 3 canopy 
types. So as a flrsv attempt at equating spectral radiance <nd1c1es and 
live aerial biomass, we combined data from each c f f he morphologic 
groups Into a single model. These data were gathered from stands of 
Typha angustlfol la , Spartlna alternlflora , Sclrpus olneyl and a broadleaf 
mixture normally dominated by AcnTda cannabina or Hibi s cus moscheutos . 

The areas were sampled monthly Trom May through August with radiance" and 
biomass data being gathered identical to the brackish marsh samples 
described in tne methods. 

Bracki s h Harsh - Biomass Estimation 

The brackish marsh sampled tor biomass estimation was diverse In 
terms of plant species and seasonally dynamic with respect to dominance. 
Table 3 lists the abundance (expressed as a percent of the total number 
of samples in which the olant species occurred) and dominance (expressed 
as a percent of the total number of samples in which the plant species 
was dominant in terms of biomass) of each plant species encountered 
during the study. The three most abundant species were S partina p atens , 
Eleocharis sp. and Acnida cannabina . Grass species dominated most 
samples during all three 'sampling periods. There was a slight Increase 
In dominance by broadleaf species In August when most broadleaf plants 
reached peak biomass. 

Simple and multiple regression models equating total live biomass 
(all species within sample) and the spectral radiance Indlcles with 
other canopy descriptors are found In Table 4. Multiple regression 
models were included for comparison because the percent of total live 
biomass which was broadleaf biomass (B) and total live biomass expressed 
as a percent of live plus dead biomass (P) were consistently implicated 
as important parameters for best fit models. Considering the very 
different spectral characteristics of broadleaf canopies relative to the 
gramineous and leafless canopies, and the potential importance of dead 
vegetation in determining spectral radiance in gramineous and leafless 
canopies, it was not surprising that these additional parameters (B,P) 
significantly improved the linear fit (r 2 values) of the models. 

Harvest and spectral radiance Index estimates of total live biomass 
for the entire growing season are compared in Table 4. In the case of 
the vegetation index, the inclusion of B and P reduced the accuracy of 
the biomass estimate, whereas with the infrared index these two parameters 
apparently improved the biomass estimate but not to any significant 
degree. The annual live biomass means predicted with radiance index 
data were very similar to the live biomass estimated by harvesting. A 
paired T-test suggested that most harvest and predicted means were not 
statistically separable at the 0.05 level. 
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Table 4 


REGRESSION MODELS EQUATING RADIANCE INDICES 
AND BRACKISH MARSH CANOPY PARAMETERS. 


Model 3 

T = (3278) VI - 1939 

T = (4319) VI - ( 7 . 1 2 ) B - 2682 

T 3 (4025) VI - (7.16)B + (1.85)P - 2562 

T * (3877) II - 2181 

T = (4242) II - (4.06)8 - 2402 

T = (3902) II - (4.41)8 + (2.35)P - 2294 


Annual 

Live Biomass (gdw 

m“ 2 ) 

Harvest^ 

Predicted 

Difference 

691(68) 

647(33) 

45 

69 ', ( 68 ) 

550 ( 49 ) 

142 * 

691(68) 

553(49) 

138* 

691(68) 

727(45) 

36 

691(68) 

680(51) 

12 

691(68) 

674(52) 

18 


a T 3 total live biomass (gdw m" 2 ), VI = vegetation index, II 3 infrared index, 

B = live biomass of broadleaf species expressed as a percent of total live biomass, 
P = total live biomass expressed as percent of total live and dead biomass, 
n for all models is 57. 

k Values are the mean and one standard error, in parenthesis, of 47 samples. 
Difference 3 the numerical difference between the harvest and predicted means. 

An asterisk indicates the difference between means was statistically significant 
at the 0.05 level. 

The regression models were then used to estimate biomass at each of 
the 3 sampling dates (Table 5). Although there was good agreement 
between annual harvest and predicted biomass means,, substantial deviation 
was noted over the growing season. Vegetation index biomass estimates 
were very good during dune and August but unacceptable during September. 

The infrared index only produced acceptable biomass estimates during 
August. The failure of the models to yield good biomass estimates 
during September was probably related to the rapid senescence of the 
marsh plants at this time. The brackish marsh, unlike the salt marsh, 
contains annual broadleaf plants which after seed production senesce 
rapidly. The larger amount of dead vegetation in the canopy and the 
senescence of broadleaf leaves altered the canopy spectral characteristics 
sufficiently to invalidate the predictive models at this time of year. 

The infrared index estimates of biomass were high in June and low in 
September. The high percentage of moisture in young, productive leaf 
tissue (June) as opposed to the low percentage of moisture in dead or 
senescing leaf tissue (September) probably contributed to the observed 
fluctuations in biomass estimation for the infrared index. 
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Table 5 

COMPARISON OF LIVE BIOMASS ESTIMATES FOR A BRACKISH MARSH BY MONTH 



DISCUSSION 


Other Investigators have found good correlations between marsh 
plant biomass and reflectance In the red and near Infrared spectral 
regions (Drake 1976, Bartlett and Klemas 1981). Me have reaffirmed this 
relationship for S. a 1 tern 1 flora using thematic mapper bands 3 and 4. 

Previous studies Tn salt marshes (Buaci and Milton 1982), pasture land 
(Curran 1982) and In agricultural vegetation (Tucker et al. 1981b) 
suggest a strong relationship between plant biomass and the combination 
of red and near infrared spectral data. These studies also imply that 
biomass could be predicted using spectral data. Very few investigators 
have actually tested models equating biomass and spectral data using a 
data set independent of the data set used to establish the model. 

Jensen (1980) working with the salt marsh shrub ( Halimione portulacoides ) 
and Curran (1980) working with pasture vegetation have done this with 
reasonable success. Our data also suggest that spectral radiance models 
are useful for nondestructive estimates of salt marsh biomass. 

The brackish marsh regression models predicted biomass well during 
the peak of the growing season. The models presented represent an 
oversimplification of the complex interactions of live and dead vegetation, 
horizontal and vertical leaf area index, and soil reflectance which are 
compressed and treated as a single reflecting surface. Brackish marsh 
canopies are normally much deeper and much more diverse in terms of 
plant morphologies thar\ salt marsh canopies. It Is, therefore, very 
encouraging when representatives of the 3 most common canopy types 
measured in relatively pure stands can be combined Into a regression 
model which yields good estimates of biomass for mixed plant stands. 

The hypothesis "that plant morphology (canopy. type) is more important ir. 
determining measured spectral radiance than consideration of each particular 
plant species, appears to be valid. For example, Tvpha canopies were 
Included as input to the predictive regression model s, yet in the data 
set used to test the model, no Typha was present. On the other hand, 
many plant species were found 1 n the test data set (notably S. patens , 
Eleocharis sp . , spicata and punctatum ) which were not Included In 

the model development data set. The dissimilarities in species composition 
between the model development and model testing data sets apparently had 
little effect upon the outcome of live biomass predictions from radiance 
Indicles. 

The vegetation index was usually a better spectral transformation 
for biomass estimation than was the Infrared Index. The vegetation 
Index has been well established as a useful transformation for monitoring 
vegetation (Tucker 1979, Tucker et al. 1979a, b) but the infrared index 
is relatively new. Kimes et al. (1981) and Markham et al. (1981) were 
among the first to work with the middle infrared band (TM5) using hand- 
held radiometry. In both studies, they concluded that the middle infrared 
band contained the same information as the red band. They did acknowledge 
that their sampling precluded any water stress and that the middle 
Infrared band should undergo additional testing when changes In leaf 
moisture were expected. More recent work by Hardlsky et al. (1983b) 
suggests that the infrared index may have moisture detection capabilities 
and therefore, will probably be useful in wetland systems for discriminating 


vegetation exposed to varying soil salinity conditions. We believe that 
additional research must be conducted to determine the worth of the 
Infrared Index for vegetation monitoring. The large seasonal variation 
In biomass predictions In the brackish marsh by the Infrared Index 
suggest moisture content or the moisture contrast between live and dead 
vegetation may have contributed to the results. If our observations are 
a result of canopy moisture differences, this would Indicate a larger 
seasonal fluctuation In canopy moisture relative to biomass than In 
chlorophyll content relative to biomass. This could potentially be a 
very useful tool in monitoring wetlands vegetation. 

CONCLUSIONS 

Regression models equating total live biomass and spectral radiance 
indicies were developed and tested for salt and brackish marsh vegetation. 
Comparisons of biomass predicted using spectral radiance indicies and 
biomass estimated by traditional harvest techniques were very similar 
for S. alterniflora biomass. The vegetation Index was slightly better 
for predicting biomass than was the infrared Index. The best agreement 
between predicted and harvested biomass occurred during August (at peak 
biomass) with a considerable amount of variability at other times of the 
year. Annual net aerial’ primary production estimates were also very 
similar using either predicted or harvest biomass estimates. 

Three morphologically distinct canopy types were Identified In the 
brackish marsh vegetation. Data gathered from broadleaf, gramineous and 
leafless canopies were combined Into a single regression model for 
estimating brackish marsh biomass. The models provided similar estimates 
of biomass compared to harvest estimates of biomass during June and 
August for the vegetation Index and only during August for the Infrared 
index. Percent broadleaf biomass and percent live biomass were Identified 
as being important parameters for determining total biomass with spectral 
data. In practice, these additional parameters did 1 little to improve 
live biomass predictions over the live biomass predictions using only 
the spectral radiance index. 

The thematic mapper bands 3, 4 and 5 used in this study successfully 
provided the spectral information necessary for nondestructive biomass 
estimates in coastal marshes. The ground-based radiometric technique 
described can provide the data necessary for estimates of productivity 
for some marsh systems. Spectral data gathered with hand-held radiometers 
from low altitude aircraft and thematic mapper simulator data are presently 
being tested using the models described In this paper. Preliminary 
results indicate that with an atmospheric correction, the models work 
well with spectral data gathered from higher altitude platforms. It 
seems plausible that the models presented here can be modified for use 
with thematic mapper spectral data. 
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